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1 Motivation 

Effects of the bottom-quark mass, mj, have been already 
noticed in the early tests El of the flavour independence 
of the strong coupling constant, a s , in e + e~-annihilation 
at the Z°-peak. They became very significant in the 
final aoalysis, which included millions of hadronic Z°- 
decays Q'B'Qta. For example, if the b-quark mass is ne- 
glected, the ratio a b s /a l g 9ht , where a b s measured from 
hadronic events with b-quarks in the final state, and 
a hgkt f rom events with light quarks (uds), is shifted from 
one El by 8%. Thus, high LEP/SLC precision requires an 
accurate account for the heavy quark mass^j in the theo- 
retical predictions for the e + e _ -annihilation into jets at 
the Z°-pole. 

The quark mass effects in the Z° decays were dis- 
cussed in the literature The leading order (LO) com- 
plete Monte-Carlo calculation for e + e~ ^ 3jets,4jets 
with massive quarks was first done in u. Later, moti- 
vated by the remarkable sensitivity of the three-jet ob- 
servables to the value of the quark mass, the possibility of 
the determination of mt at LEP, assuming universality of 
strong interactions, was consideredtru. This question was 
analyzed in detail intl, where the necessity of the next- 
to-leading order (NLO) calculation for the measurements 
of the mb was also emphasized. 

The NLO calculations for the process e + e~ — > 3jets, 
with complete effects of the quatjk mass, were performed 
independently by three groupsEjIiilO. These predictions 
are in agreement with each other and were successfully 
used in th&jneasurements of the b-quark mass far above 
thresholdaEl and in the precision tests of the universality 

a It is mainly related to the bottom-quark. Effects of the charm- 
quark mass are smaller, roughly by the factor mj?/mj\ 



of the strong interactions iHHi at the Z°-pole^. 

In this talk we make a short review of next-to-leading 
order predictions for e + e~ — ■> 3jets including effects of 
the quark mass and its applications at the Z°-peal|]. 

2 Why are the b-quark mass effects are signifi- 
cant in Z° — > 3jets? 

It might seem surprising that at the Z°-pole, where the 
relevant scale is the Z°-boson mass, effects of the quark 
masses can not be neglected. Indeed, if one considers 
inclusive width of the decay Z° — > bb, according to the 
Kinoshita-Lee-Nauenberg theorem there are no mass sin- 
gularities, and the only way m& enters calculations^ is in 
the ratio W^(Mz)jM\. Therefore, quark mass effects in 
the best-measured observable for b-quarks are negligibly 
small ~ fO -3 . 

But the situation is different in more exclusive pro- 
cesses. Let's consider the decay Z° — > bbg, which con- 
tributes to the three-jet final state at the LO. This pro- 
cess has an infrared singularity in the limit when the en- 
ergy of radiated gluon energy aproaches zero. To make a 
physical prediction one has to introduce kinematical cuts. 
In the e + e~-annihilation this is usually done bjuapply- 
ing one of the so-called jet clustering algorithms Ej. The 
phase-space for Z° — ► bbg is split into two parts, two-jet 
and three-jet one, and this separation is defined by the 
jet-resolution parameter, y c . Therefore, instead of two 

b Due to the large correlation between a s and mj, either rnj, or 
a s was treated as a free parameter. The value of another one was 
taken from other measurements. 

c We would like to note that elements of these calculations can be 
also applied for the e+e — — * tt + • • •. 

d Taking the MS running quark mass one includes ifchpi principal 
part of the NLO QCD corrections to the total widthEjij. 



scales in the inclusive process, Mz and mb, one has here 
in addition the new scale, ^fy~ c Mz- The jet-resolution 
parameter can be rather small, in the range 10~ 2 — 10~ 3 
and effects of the quark mass appears as m 2 / \^fy~ c M z) 2 
and can reach several percents. 



3 The three-jet observable 

The convenient observable for studies of the mass effects 
in the three-jet final state is defined as follows 



r>bd 1 3] 

H 3 = - 



(i) 



= 1 + n (b (y c ,r b ) + ^&i(y c ,r 6 )) + 0(a 2 ) 

where r|- and T q are the three-jet and the total de- 
cay widths of the Z°-boson into a quark pair of flavour 
q, ?'b = m 2 /M z . Both the LO function, bo, and the 
NLO function, b\, depend on the jet-clustering algorithm. 
Note that although the leading ^-dependence is factor- 
ized for convenience, the above expression is not an ex- 
pansion in rv 

This observable has both experimental and theoret- 
ical advantages. It is a relative quantity, therefore many 
experimental uncertainties due to the normalization drop 
out. In addition, in the ratios T^/T q the bulk of elec- 
troweak corrections is cancelled. The double ratio ([I]) 
was measured by DELPHI □ and used for the determina- 
tion of the rtib- In practice one uses normalization with 
respect to all light flavours, uds. Such quantity then 
differs from R 3 d , mainly due to the contribution of the 
triangle diagrams t3. The difference is, however, very 
small numerically Similar ratio normalized on udsc was 
considered mil. 

Another observable is the differential two-jet rate, 
L>2 defined as 



D bd = 



[T b 2j (y c + Ajfe/2) - r b 2j (y c - Ay c /2)}/T b 
[T%{y c + Ay c /2) - T%{y c - Ay c /2)]/T d 



(2) 



with Ay c taken to be sufficiently small. The two-jet 
width is obtained from the relation: r|„- = T q 



where T\ - is the four-jet width. 



T q - 



-V q 
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nt event shapes observables were consid- 
in the literature as well. 



4 The leading order calculation 

The LO contribution to the decay Z° — > 3jets is given 
by the process Z° — > bbg. At the LO we can not spec- 
ify what value of mb should be taken in the calcula- 
tions: all definitions of the quark mass are equivalent 
(the difference is due to the higher orders in a s ). One 



can usep| for example, the pole mass Mb w 4.6GeV or 
the AfS'-running mass rfib{^) at any scale relevant to the 
problem, mb < /i < Mz, with mb(mb) w A.lSGeV and 
Mb(M z ) w 2.83GeV. 
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Figure 1: The ratio R^ d (Jade) as a function of y c 
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Figure 2: The ratio R^ d (Durham) as a function of y c . 

The xesults for the ratio R b 3 d , Eq.(0), for Jade0 and 
Durham t3 schemes are presented by solid curves in the 
Figs.0J|. The upper curve corresponds to mb = 2.83GeV, 
the lower one to nib = 4:.6GeV. The difference between 
the two curves gives the size of the uncertainty of the LO 
prediction. To improve the situation the NLO calculation 
is necessary. 

It is worth mentioning that the residual mass de- 
pendence in the function bo(y c , rb) is small in the region 
0.005 < y c < 0.1, especially for the Durham scheme. The 



simple interpolatio: 
be used in practice! £|, ' J 



- &1 2 ^ In 2 y c can 



e The values of the b-quark mass«s are taken from the recent sum 
rules and lattice QCD analyses Ej of the T and B mesons spectra. 



5 The next-to-leading order corrections 

At the NLO there are two different contributions. One 
comes from one-loop corrections to the three-parton de- 
cay, Z° — ► bbg. Another one comes from tree- level four- 
parton decays Z° — > bbgg and Z° — > bbqq, q = u,d, s, c, b 
integrated over the three-jet region of the four-parton 
phase-space. The main difficulty of the NLO calcula- 
tion is the presence of ultraviolet (UV) and infrared (IR) 
divergences at the intermediate stages. The UV diver- 
gences in the one-loop part are removed by renormaliza- 
tion of the QCD parameters. The IR singularities in the 
virtual part are due to massless gluons in the loops. In 
the four-parton process they appear when one of the glu- 
ons is soft or two gluons are collinear. The sum of the 
virtual and tree-level contributions is, however, IR finite. 
In addition, the finite quark mass makes the NLO calcu- 
lation technically much more involved, pComparing to the 
massless one, known since many years ErEU. 

The one-loop part is calculated analytically and di- 
mensional regularization is used to regularize both UV 
and IR divergences. The IR singularities appear as sim- 
ple and double poles ha e, D = 4 — 2e, where D is space- 
time dimension. The singular part is proportional to the 
tree-level transition probability and it is cancelled by the 
IR-singular part of the four-parton contribution. 

There are severaLmathods of analytical cancelation 
of IR singularities. InEJEJ the so-called slicing methodEj 
was used. In this case the analytical integration over 
a thin slice at the border of phase-space is performed in 
D-dimensions. The integration over the rest of the phase- 
space, defined by a particular jet-algorithm, is done—nu- 
merically in D = 4. In the third NLO calculation E3, a 
different approach, the so-called subtraction method (see 
e.g. E3a) was used. 

We would like to stress that the structure of the 
NLO result for the Z° — > 3jets in the massive case is 
completely different from the massless one EH. In the 
massive case the collinear divergences associated with the 
gluon radiation from the quarks, are softened into lnr& 
and only collinear divergences due to gluon-gluon split- 
ting remain. Therefore, the test performed inE£l by recov- 
ering the massless limit from the result with finite mass, 
is a rather non-trivial one. 

In contrast to the LO function, bo, the NLO frujpj 
tion b\ has a significant residual mass dependence 
and the phenomenological interpolation, e.g. for Durham 
scheme, can be chosen as follows 



do) 



b^' lny c + bY' lnr^ . 



(2) 



(3) 



The approximation Eq.([|) works well in the region 0.01 < 
y c < 0.1 and extra powers of In r\> and/or mixed terms 
In ?'b In y c do not improve its quality. 



In the NLO calculations one can, and have to, spec- 
ify the quark mass definition. It turned out that techni- 
cally it is simpler to use a mixed renormalization scheme 
with on-shell definition for the quark mass and MS def- 
inition for the strong coupling. In this case observables 
are originally expressed in terms of the pole mass. This 
definition of the quark mass can be perfectly used in per- 
turbation theory. However, in contrast to the pole mass 
in QED, the quark pole mass is not a physical parameter. 
The non-perturbative corrections to the quark self-energy 
bring an unavoidable ambiguity of order « 300Af eV 
(hadron size) to the physical position of the pole of the 
quark propagatorcil Above the quark production thresh- 
old, it is natural to use the running mass definition (we 
use MS). The advantage of this definition is that ra^/j) 
can be used for any scale, jU » m&. The pole, Mb, and the 
running masses of the quark are perturbatively related 



M b = mb(p) 



1 



Ots_ 

n 



In- 



(4) 



Although this relation is known to higher orders in a s , we 
use its one-loop version to pass from the pole mass to the 
running one, which is consistent with NLO calculations. 
Substituting Eq.(|^) into Eq.(Q) we have 



Rl a (y c ,m b (n),Li) = 



1+F 6 (m) 



a s (n) 



M 2 

- bl -2b \n^ 



(5) 



with b 1 = bi + b (8/3 - 21nr 6 ) and r b = ml/M%. 

To relate quark masses at different scales, we use the 
one-loop renormalization group equation, e.g. rrib(p) = 



m b (M z ) 
with rif 



a s (Mz) 

5. 



so 



where 70 = 2 and (3q = 11 



6 Numerical results and discussion. 

In Figs|l],|2| we show NLO results, obtained by using 
Eq.(|^) parameterized in terms of rrlb(Mz) = 2.83GeV for 
different values of the renormalization scale /i = lOGeU 

- dashed line, fi = SOGeV - dotted line and /i = Mz 

- dashed-dotted line. As expected, the NLO curve for 
large scale is closer to LO curve for rnb(Mz), while NLO 
result for smaller scale is_closer to the LO one for Mb- 
Like in the massless caseO, the NLO corrections for the 
Jade scheme are larger than in the Durham and depend 
significantly on y c for small values, y c < 0.025. 

The /i-dependence is a reflection of the fixed order 
perturbative calculation and the size of the variation 
of the observable is usually used to estimate higher or- 
der corrections. The study of the /x-dependence for the 
Durham scheme is presented in Fig.||. The R^ d is shown 
as a function of scale p, for the fixed value of y c = 0.02. 



The lower, solid curve is obtained by using Eq. (|l|) , i.e. pa- 
rameterized in terms of the pole quark mass. In this case 
the scale dependence is due to the renormalized coupling 
constant, a s . Other curves show /x-dependence when R^ d 
is parameterized in terms of the running mass, rnb(Mz), 
Eq. (||) , with different definitions of the quark mass used 
in the logarithms: M&, rrib(Mz), mf,(/x). 
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Figure 3: The ratio R^ d (Durham) as a function of /x. 

The conservative estimate^ of the theoretical error for 
the i?3 d is to take the whole spread given by the curves 
in Figj^. The uncertainty in R^ d induces an error in 
m b {M z ): AR^ d = 0.004 -► Am b ~ 0.23GeV, which is, 
however, below current experimental errors, dominated 
by fragmentation. 

In the FigsJ^J^ we show results for the fi&ed value of 
scale, jj, = Mz , for the Cambridge algorithmic, a recently 
proposed modification of the Durham scheme. 
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Figure 4: The ratio R^ d (Cambridge) as a function of y c - 

It is remarkable that in this scheme and for two 
observables, R^ d and two-jet differential ratio, Di (see 
Eq.(|])), the LO result for the running mass, mb(Mz) is 
very close to the NLO one. In a sense, it is similar to 
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the total width: the main radiative corrections are taken 
into account by the running of the QCD parameters to 
the M^-scale. Note also that for small, but still reason- 
able value of y c sa 0.01, the mass effects in D2 are as 
large as 10% (although one has to remember that this 
is a differential rate and statistical errors are larger here 
too). 
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Figure 5: The ratio (Cambridge) as a function of y c . 

The scale dependence for the Cambridge scheme is 
shown in Fig.^J for fixed y c = 0.01. The result is also 
quite remarkable. The ratio expressed in terms of run- 
ning mass is very stable with respect to scale variation. 
This behavior in the Cambridge scheme is rather promis= 
ing with respect to improvements of the DELPHI resultB 
for mb(Mz)- We refer tocB for detailed discussion of the 
NLO predictions in this new scheme and to E3 for the 
first analysis of the LEP data applying the Cambridge 
algorithm. 
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f The Jade scheme gives significantly larger error, see FigJI]. 



Figure 6: The ratio R^ d (Cambridge) as a function of fi. 



7 Conclusions 

Three-jets observables at the Z°-peak (jet rates, differ- 
ential jet-rates, event-shape observables etc.) have sig- 



nificant mass effects ranging up to the 10% depending 
on the observable, jet algorithm and the value of the jet 
resolution parameter. This requires accurate theoretical 
input, including mass effects, for tests of the flavour in- 
dependence of the strong interactions and measurements 
of the bottom-quark mass. 

In the last years an important progress was done in 
the description of the decay Z° into three-jet with mas- 
sive quarks. The next-to- leading calculations have been 
done by three groups and have been successfully used in 
the analysis of the LEP and SLC data. The NLO correc- 
tions are in the range of 1 — 3% and are within the ex- 
perimental reach. Further studies of different observables 
and different jet-algorithms are oriented on the reduction 
of the theoretical uncertainty. One good candidate might 
be the Cambridge jet- algorithm, where the NLO correc- 
tions are particularly small and where the predictions in 
terms of the running mass, mb(Mz) are particularly sta- 
ble with respect to the variation of the renormalization 
scale. 
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